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ABSTRACT 

We present mid- infrared observations of 18 sources from a sample of 21 2; ~ 2 
radio-intermediate obscured quasars. The mid-infrared spectra of the sources 
are continuum dominated, and 12 sources show deep sihcate absorption with 
r9,7~l-2. Combining mid-infrared and optical spectra, we achieve 86% spec- 
troscopic completeness which allows us to confirm that most (68^22%) -2^2 
radio-intermediate quasars are obscured. The new spectra also prove that many 
high-redshift type-2 quasars do not show any rest-frame ultraviolet emission lines. 
From the 18 individual mid-infrared spectra, we classify most of the sources into 
three subsamples: those with hints of the 7.7 and 6.2 //m polyaromatic hydrocar- 
bons (3/18 sources show PAHs, subsample A), those with an excess of emission 
around 8 /im but no hint of the 6.2 /im PAH (7/18 cases, subsample B) and 
pure-continuum sources with no visible excess (4/18 sources, subsample C). The 
remaining 4/18 sources have spectra that are featureless or too noisy for any fea- 
tures to be visible. In subsample A, averaging the spectra leads to a statistical 
detection of both 6.2 and 7.7 /im PAHs over the continuum, with the strength of 
the 7.7 /im PAH comparable to that of submillimetre-selected galaxies (SMGs) at 
similar redshifts. These sources are in a phase of coeval growth of a supermassive 
black hole and a host galaxy. 

Subject headings: galaxies:nuclei - galaxies:active - infrared: galaxies - quasars:general 
- galaxies:starbursts - galaxies:high-redshift 



1. Introduction 

The correlations between galaxy and central supermassive black h ole properties found 



in the local Universe (IMagorrian et al.lll998l : iFerrarese fc MerrittI 120001 ) have led to notions 
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that the growth of both must be intimately related. A fraction of the gas fuelling star 
formation is b elieved to rea ch the black hole, which will grow by accretion during an obscured 
quasar phase (lFabianlll999l ). and is expected to develop winds or jets that eventually becona e 
powerful enough to expel the gas, preventing any further star formation (jSilk &: Reeslll998l ). 
This feedback proces s is needed to avoi d forming more stars than are seen around local 
dormant black holes (IBenson et al.l 120031 ) . but few compelling examples of this phase have 
been reported to date. 

Quasars are believed to be viewed as obscured or unobscured depending on the ori- 
entation of the line of sight wit h respect to the axis of symmetry and the dusty "torus" 
surrounding the accretion disk (lAntonuccil Il993l ) . yet the population of powerful obscured 
quasars had remained elusive. Since most of the growth of supermassive black holes is 
obscured by dust, optical surveys miss this phase completely, and X-ray surveys miss a sig- 
nificant fraction due to associated gas column. Recently, mid-infrared surveys have been able 
to find obscured active galactic nuclei (AGN) in large numbers, but care must be taken to 
separate AGN from star-forming gala xies. This can be achieved by making use of t he particu- 



lar mid-infrared colours of AGN (e.g. iLacy et al. 



comb ination of mid-infrared and radio data (e.g. 



2004 



Alonso-Herrero et al. 



Martmez-Sansigre et al 



200 



i 



20061) or using a 



Donley et al. 



20051 ). Indeed, from thes e surveys it seems t hat a large fraction of obsc u red AGN are un- 



detected in X-ravs fe.g. iDonlev et al.l l2005l : I Alonso-Herrero et al.l l2006l : iLacy et al.l 12007 



Martmez-Sansigre et al.ll2007l ). and that mid- infrared selected samples give a more complete 
picture than X-ray selected samples. Some obscured quasars, however, do not show even 
the narro w emission lines expected if th e centra l regions are obscured at certain angles by a 



torus (see 



Martmez-Sansigre et al 



narrow-line region (IBrand et al. 



200 



20051 . l2006al ). some have evidence for dust in front of the 
and radio observatio ns suggest that in some cases the 



radio jet is closely aligned with the observer's line of sight (iMartfnez-Sansigre et al.ll2006bl : 
Sajina et al.ll2007l . and Klockner et al., in prep.), observations which are incompatible with 
the scenario of all obscuration being caused by the torus around the accretion disk. 

The existence of objects which are not easily explained in the torus model has led to the 
suggestion that, in many cases, dust on scales of kpc is responsible for the obscuration. Since 
young stars are born in dusty regions, it is possible that in these sources the host galaxy is 
undergoing an episode of star formation. Mid-infrared spectroscopy gives a clean method of 
looking for star formation in the form of emission from complex molecules like polycyclic aro- 
matic hydrocarbons (PAHs, e.g. lAitken fc Rochelll985l : iRoche et al.lll99ll ). These molecules 
are excited by ultraviolet photons from young stars, but destroyed by the hard rad iation of 
an AGN out to ~ 100 pc - 1 kpc scales (e.g. IVoitlll99ll : ISiebenmorgen et al.l 12004 ). Thus, 
detection of the PAH emission bands can only be attributed to star-formation and not to 
excitation by an AGN. 
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If obscured and unobscured quasars are identical types of objects viewed at a different 
orientation, then tlie properties of their host galaxies, such as the PAH lum inosity, should be 
identical on average. In a detailed study of local AGN, IClavel et al.l (I2OOOI ) found that indeed 
the PAH luminosities of Seyfert 1 and Seyfert 2 AGN are statistically indistinguishable, and 
that the observed differences between the Seyfert 1 and 2 spectra are due to the different 
contrast between the PAHs and the continuum, the latter being weaker for Seyfert 2s. 

In this letter we present the results of mid-infrared spectroscopy of a sample of radio- 
intermediate obscured (type-2) quasars, which show either narrow lines or no emission lines 
in their optical (rest-frame ultraviolet) spectra. Throughout this paper we assume a ACDM 
cosmology with h = i/o/(100 km s'^ Mpc~^) = 0.7; Qm = 0.3; Qa = 0.7. 



Observations and data reduction 



Our sample of radio-intermediate obscured quasars, first presented in lMartmez-Sansigre et al 



(l2005l ). was observed using the Infrare d Spectrograph (IR S) instrument (IHouck et al.ll2004l ) 



on board the Spitzer Space Telescope (IWerner et al.ll2004l . GO-30634 and archival data). At 



^ ^ 1, the brightest PAH lines as well as the silicate features are all redshifted out of the 
range covered by the Short-Low module (5.2-14.5 fim). Therefore, only the low resolution, 
long wavelength detector was used (Long-Low), with both orders, giving coverage in the 
spectral range 14-38 fim. However, due to the high noise above 35 fim, we only make use of 
the range 14-35 /im. The number of cycles each source was observed for is shown in Table [1] 
Each cycle involves two nodded positions, with an individual cycle duration of 120 s at each 
position. Note that three sources were not observed due to their high-redshifts (AMS03, 
AMS12 and AMS16). 

The basic calibration data (BCD) from the Spitzer Science Center IRS pipeline were 
then combined by subtracting nodded pairs, and median combining the pairs. The number 
of nodded pairs combined in this way varies between 2 and 10, depending on the source (see 
Table 1). Both mean and median combining were tried, and for most sources the resulting 
spectra were almost indistinguishable. Mean combining lead to slightly better signal-to-noise 
ratio (SNR) spectra in the smooth regions, but median combining generally proved better 
at removing spikes of noise. Given the low SNR for most of our spectra, median combining 
was chosen. 

The 2-dimensional spectra were cleaned of rogue pixels using irscleanEI, and 1- 



http: / /ssc. spitzer. caltech.edu/postbcd/softwarehome. html 
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dimensional spectra were extracted using the optimal extraction in SPICE^. For the optimal 
extraction, an uncertainty file was first created from the frame to frame variance of each 
pixel, using the 440 BCD files of our own GO program. This uncertainty file was scaled to 
the correct exposure time for each source, and is overlayed in Figure 1. All 18 sources are 
detected in the reduced data. 



3. Mid-infrared spectra 

3.1. Spectral characteristics 

The mid-infrared spectra allow us to obtain spectroscopic redshifts for a further 8 
sources, confirming that 17 out of 21 sources in our sample lie at z > 1.6. Contami- 
nation by pure starbursts is unlikely, since the mid-infrared spectra are characteristic of 
AGN-dominated sources. 

Figure 1 shows the 18 individual spectra. Silicate absorption can be found in most 
of the sources. In two cases (AMS05, see Appendix Rl AMS17), the lack of this feature 
can be explained due to their high redshift, so the feature is expected to be centred at and 
observed wavelength Aobs ^ 37 fim. In other sources the low SNR means it is difficult to 
see clearly this feature, although hints of it are seen in the 2D and ID spectra. In addition, 
from Figure 1, it can be seen that AMSOl could be at z ~ 2.1 and have a shallow silicate 
feature, but we do not consider this redshift secure so we do not quote it in Table 1. 

In total, 12/21 sources show no emission lines in the optical (rest-frame ultraviolet, see 
also Appendix A). The mid- infrared spectra of these 12 sources are continuum-dominated 
and 8 sources show the silicate absorption feature. These mid-infrared properties are char- 
acteristic of AGN. In fact 6 sources showing no rest-frame ultraviolet lines are at redshifts 
where the Ly a falls in the optical range (that is z > 1.7), and therefore have no Ly a 
detection down to Si^ya <2 x 10~^° W m~^. We estimate the absorption to the broad-line 
region, at 1216 A to be Ai^ya ~ 60 — 120 magnitudes (see Section 13.21) . The extinction 
towards the narrow-line region is expected to be significantly lower. However, the Ly a 
line suffers from resonant scattering, so modest amounts of dust mixed with the narrow-line 
region could still severely extinguish the lin e emission. Indeed, the Ly a profile of AMS03 



Figure 6 of iMartmez-Sansigre et al.ll2006a|l. is ve r y sini ilar to the simulated spectrum of a 



young dusty galaxy by lLaursen fc Sommer-LarsenI (120071 . Figure 2), which includes resonant 



scattering. AMS03 is close to the limit of our sensitivity, so sources where the narrow-line 
region is dustier or where the Ly a line is intrinsically less luminous than in AMS03, will 
not show Ly a in our optical spectra. 
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We find 10 objects with an excess emission around rest-frame 8 /im, wliicli at first sug- 
gests the 7.7 /im PAH emission band. These objects are AMS02, AMS04, AMS06, AMS08, 
AMS09, AMSll, AMS15, AMS17,AMS18 and AMS21. Care must be taken, however, since 
for heavily absorbed sources, there is an apparent excess at 8 /xm. This is due to the mini- 
mum in optical depth between the absorption features due to sihcate, water ice (centred at 
6.0 Atm) and hydroge nated amorphous carbons (at 6.9 and 7.3 /xm, for a detailed discussion 



see 



Spoon et al.ll2002l ). At first sight, this excess emission around 8 /xm can mimic the 7.7 /im 
PAH, but it is wider than the 7.7 /im PAH, and it is not accompanied by the 6.2 /im PAH. 
In 3 sources (AMSll, AMS17 and AMS21), there is a hint of another emission band corre- 
sponding to the 6.2 /im PAH. A mean spectrum of these 3 sources (hereafter subsample A) 
is shown in Figure 2(A), where the 6.2 and 7.7 /im PAHs are detected. The mean spectrum 
of the remaining 7 sources with a maximum around 8 /im, but no clear hint of the 6.2 /im 
PAH (subsample B), is shown in Figure 2(B). Figure 2(C) shows the mean spectrum of sub- 
sample C, consisting of 4 sources (AMS05, AMS13, AMS14, AMS19) continuum dominated 
with clear silicate absorption, but no hint of excess emission around 7.7-8 /im. For reference. 
Figure 2 has the mid-infrared spectra of severa l low-redshift ultra luminous infrared galaxies 



overlayed (ULIRGs, with Ljr > 10^^ Lq, from lArmus et all 120071 : ISpoon et al.ll2007r i. These 



are: IRAS 12514+1027 (hereafter IRAS 12514), Mrk 231, Mrk 273 and Arp 220. The mid- 
infrared spectrum of IRAS 12514 represents a pure-AGN, Mrk 231 and Mrk 273 represent 
AGN-starburst composites, while in the mid-infrared, Arp 220 is starburst-dominated. 

For each subsample, the mid-infrared spectrum of IRAS 12514 was used as a continuum 
template. It was normalised using the mean luminosity density of the 6.5-7.5 /im region, 
and then subtract ed. The residuals are plotted in Figure 3, with the spectrum of M 82 



(jSturm et al.ll2000l ) overlayed as a reference for the location of the PAHs. Table [2] summarises 



the mean PAH luminosities or limits for each subsample. 

In subsample A, there is a detection of the 6.2 /xm PAH (3.4 a, see Table ^ together 
with the 7.7 /im PAH (6.4 a). The noise is calculated in the 5.0-8.5 /im region, (note 
this is lower than in the 8.5-11.5 /im region), and as is mentioned in Table [21 it is the 
noise in the 5.0-8.5 /im region that is used to estimate the significance of the detections. 
The uncertainty in the line flux additionally includes the uncertainty in the continuum 
subtraction. In subsample A the presence of PAHs is an indicator for ongoing starformation 
in the host galaxy. 

Subsample B shows an excess around 8 /im but no hint of the 6.2 /im line. The shape of 
the excess, and the lack of the 6.2 /im PAH, strongly suggest this excess at 8 /im is continuum 
and not the 7.7 /im PAH. Thus, in subsample B we conclude there is no evidence for PAHs, 
within our sensitivity, and that the maximum at 8 /im is probably continuum. Subsample C 



- 6 - 



shows no hint of any PAH or any excess at 8 fim, and indeed the residuals in Figure 3 (C) are 
consistent with noise. The mid-infrared spectrum of subsample C is, within our sensitivity, 
indistinguishable from that of IRAS 12514. 



Wilman et al.l (120031 ) obtained an X-ray detection of IRAS 12514, and estimated the 
total bolometric luminosity to be Lboi ~1.6xl0^^ Lq, with a Compton-thick quasar and a 
powerful starburst contributing comparable fractions to Lboi (a source is Compton-thick if 
the absorbing column density, iVn, is > l/cxr =1.5 x 10^^ m~^, where ax is the Thomson 



electron scattering cross-section). The mid-infrared spectrum of IRAS 12514, however, is 
A GN-dominated, wit h a very low value of the 6.2 fim PAH equivalent width (see Figure 1 
of ISpoon et al.l 120071 ). Thus, although the mid-infrared spectra of samples B and C show 
no signs of starburst activity, vigorous ongoing star-formation in the host galaxy cannot be 
ruled out. 



3.2. Derived extinctions 



The spectra are noisy in the silicate region, and from Figure 3 it is clear that the derived 
depth of the absorption feature, tqj, will vary vastly depending on the exact wavelength at 



which it is measured and the choice of anchor for the continuum (see, e.g., ISpoon et al. 



20071 ). However, we can estimate rg y by comparison with the local U LIRGs. Arp 220 ha s 
rg.r = 3.3, Mrk 231 rg.r = 0.8, a nd Mrk 273 rg 7 = 1.8 (all three from lArmus et aDl2007[ l. 
while IRAS 12514 has rg.7 = 1.5 (ISpoon et al.ll2007l ). For subsample A, in the silicate region, 
the spectra of Mrk 273 and Mrk 231 represent lower and upper envelopes, respectively, so 
0.8 < Tg 7 < 1.8. For subsample B, IRAS 12514 provides an upper envelope and Arp 220 a 
lower one (so 1.5 < rg.7 < 3.3), while within the noise, the spectrum of Mrk 273 is a good 
fit, so rg.7 ~ 1.8. The spectrum of subsample C is very close to that of IRAS 12514, so 
rg.7 ~ 1.5. 



We can see that rg.7 is typically ~l-2. Assuming Ay = 18.5 x rg.7 (lDrainell2003l ). we 
therefore expect values of Ay in the range ~18.5-37. Using the Milky- Way dust model of 
Peil (I1992I ). one expects the extinction at 1216 A to be Ai^ya ~ 60 — 120 magnitudes. This 
is an estimate of ^Lya towards the central region (i.e., to the broad Ly a line). 



Following iMartmez-Sansigre et al.l (l2006al ) , we assume a gas-to-dust ratio of iVn = 5 x 

We compare 



2 X 10^8 m-2. 



]^g26 2 xAy, which suggests column densities ~ 9 x 10^^ 
this to the absorbing column densities of the AGN in IRAS 125 14, Mrk 231 and Mrk 273 



Balestra et all (120051 ) find iVn = 7 x 10^^ in Mrk 273, while I Wilman et al.l (120031 ) find 



Nn > 1.5 X 10^8 m-2 in IRAS 12514. For Mrk 231. lBraito et al.l (l2004f ) estimate A^^h ~ 2 x 10 



l28 
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m ^. Thus, the low-redshift ULIRGs whose mid-infrared spectra most resemble our sample 
are all AGN heavily absorbed in X-rays, or even Compton-thick. 

Indeed, the detection of PAHs in subsample A suggests the star-forming regio ns emittiri g 
the PAHs must be shielded from the X-rays, or else the PAHs would be destroyed (IVoitlll992l ). 
The line of sight of the PAHs to the central source of X-rays is not necessarily the same as 
our line of sight, but there is a strong suggestion here that many of the sou r ces iii this sample 
will be Compton-thick, as was suggested also by iMartmez-Sansigre et al.l (120071 1 . 



The mid-infrared spectra can also be used to estimate the mean bolometric luminosity of 
each subsample. Assuming the 6.5-7.5 fim region is dominated by the AGN continuum, then 
Lboi ~ 10xAL a~ 10^^ Lf7) (assu ming and unobscured quasar spectral energy distribution. 



SED, following IEI vis et al.l 1 1994 ). These sources are not, of course, unobscured quasars, but 
the effect of dust at 7 /im is small: we have estimated the range of Ay to be ~18.5-37. 
Using the Milky- Way type dust model of iPeil (119921 ). the expected transmission at 7 /im is 
~60-75% for this range of Ay, while it is ~50-55% for Small-Magellanic Cloud type dust. 
Hence, the bolometric luminosity is underestimated by at most 50%, which is smaller than 
or comparable to the intri nsic uncertainty in us ing the Elvis et al. (1994) conversion (which 



is a factor of ~ 2, see also Richards et al. 2006). Note that values of Lhoi^ 10 Lo^ are also 



estimated for the sample of iMartmez-Sansigre et al.l (120071 ). which was selected using very 
similar criteria as well as broad-band SED fitting. 



3.3. Possible selection bias? 

A source of concern is whether there is a selection effect in our sample: by being selected 
at 24 /im, it could be biased in favour of 2; ~ 2 sources with strong PAHs. To a certain extent 
this is inevitable, but we show in this section how the effect is minimal. In addition, a look 
at the spectra in Figure 1 shows such a selection effect has not heavily affected us. 

At z = 2, the 7.7 /im PAH is almost perfectly centred on the MIPS 24 fim band. We 
estimate the mean contribution to the 24 /im flux density of a L7.7 = 5 x 10^° L0 PAH (from 
Table [2]) at z = 2 to be ~100 /iJy. This can clearly be important for a sample whose flux 
density limit is 300 /iJy. Thus, at the lower flux density end, our sample could be biased in 
favour of sources with PAHs, yet we expect at most ~ 30% contribution of the PAH to the 
24 /im flux density, and for sources with S24, > 650 /iJy, the contribution of the PAH will be 
< 15% and therefore smaller than or comparable to the uncertainty in the flux density. 

Amongst the faint sources, we do see cases where the sources have been selected due to 
the maximum around 8 /im falling in the 24 /im band (e.g. AMS02, AMS15, see Figure 1). 
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Yet, as we have discussed in Section \3.1\ without detection of the 6.2 fim line, we beheve 
this 8 /im excess to be continuum and not the 7.7 fim PAH. We have only identified PAHs 
in AMSll, AMS17 and AMS21. 

In the case of AMSll, the peak of the 7.7 fim PAH is at 20 /im, but the MIPS-24 fim 
band has a spectra l response curve with a full- width half maximum between 20.8 and 26.1 fim 



(IRieke et al.l 12004 ). Hence, at most half of the 7.7 /im PAH fiux falls in the band, so the 
PAH contribution is expected to be ~ 50 /xJy. AMSll has a 24 fj,m fiux density, 5*24 = 442 



//Jy (24 /xm fiux densities from iMartfnez-Sansigre et al.l l2005l ). so subtracting the PAH 
contribution one expects S24 ~ 392 /iJy, enough to be inside our sample anyway. AMS17 
is at 2; = 3.137, so only the 6.2 /im PAH falls in the 24 /im band. It has S24 = 1134 /iJy: 
clearly the PAHs have no effect in the inclusion of this source in our sample. AMS21 is at 
z = 1.8 and the PAH falls entirely within the 24 /im band, but with 6*24 = 720 /iJy the 
PAH-subtracted flux density is 5*24 ~ 620 fiJj, enough to be included. A look at AMS21 in 
Figure 1 suggests the continuum to be slightly lower than this estimate, around 5*24 ~ 400 
/iJy, but still bright enough to be included in the sample. 

Overall, the spectra shown in Figures 1 and 2 show our sample is not contaminated by 
pure starbursts, and in any case the radio properties suggest all our sources contain AGN: 
the sources have radio luminosities, Li .4 ^ lO^^^ W Hz ~^ sr~^, several have flat- or gigahertz- 
peaked radio spectra (IMartmez-Sansigre et al.ll2006bl ). and many have radio cores detected 
(Klockner et al., in prep.). 



4. The quasar fraction revisited 



We use our new mid-infrared spectra (and t he new optical spectra i n App endix lAl) 
to revisit the quasar fraction at z ~ 2, following iMartmez-Sansigre et al.l (120051 ). Above 
2 = 1.7, Ly a is potentially visible in optical spectroscopy, while in our mid-infrared spectra, 
the 9.7 /im silicate absorption feature is visible in the approximate range 1.4 ^ z < 2.6. Thus, 
our spectroscopic completeness to both "host-obscured" and "torus-obscured" quasars should 
now be close to 100% in this range. Above z ~ 2.6 we are no longer complete to sources 
showing no emission lines in the optical spectra. Thus, we estimate the quasar fraction at 
1.7 < z < 2.6, where we expect to have close to 100% spectroscopic completeness. We note 
that objects with no narrow lines in the optical spectrum and with shallow silicate absorption 
features could still be missing from our census (e.g., AMSOl might be at z ~ 2.1 as discussed 
in Section 13711) . 



We begin by estimating the number of type-1 quasars following our 24 /im and radio 
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criteria. This is estimated by using the IWolf et al.l (120031 ) luminosity function (assuming pure 



luminosity evoluti on, PLE), together wit h a typical type-1 SED (essentially fiat in vLy in 



the infrared, from iRowan-RobinsonI Il995l ) are used to estimate the number of 1.7 < z < 2.6 



type-1 quasars meeting our 24 /im cut. The optical-to-radio correlation of ICirasuolo et al. 



(120031 ) is then used to see what fraction of these type-Is would also meet the radio cuts. We 
estimate the errors by changing from PLE, to pure density evolution (PDE). We also vary 
the mid-infrared spectral index (where oc z/""'^!^^) of the typical type-1 SED, assumed to 
be ttMiR = 1, in the range 0.87 < aMiR < 1.13, and the radio spectral index, assumed to be 
ttrad = 0.8, in the range 0.5 < c^rad < 1- The errors are then added in quadrature, with the 
difference between PDE and PLE considered as a ±lo" error, while the range spanned by 
the spectral indices is considered to be ±2o". In an area of 3.8 deg^, we predict 6.0l^ 4 such 
type-1 quasars. 

Our sample includes 4 objects with spectroscopic redshifts 1.7 < z < 2.6 from rest- 
frame ultraviolet lines in the optical spectra. A further 6 objects have 1.7 < z < 2.6 from 
the mid-infrared spectra (see Table 1). Figure 4 shows the normalised posterior probability 
distribution of the quasar fraction, g, given our data, given the modelled number of type-1 
quasars, and assuming a prior probability distribution for the quasar fraction fiat in the entire 
possible parameter range of < g < 1 (again, following Martinez-Sansigre et al. 2005). The 
solid blue line shows the quasar fraction when only the type-2 quasars with optical redshifts 
are used (4 type-2s for a predicted 6.0 type-Is): it has a modal value of 0.60 and 68% of the 
area lies in the region with q > 0.51 (note that, for this curve, the errors are so large that 
the probability of g > 1.0 is greater than zero). The dashed red line shows the posterior 
distribution for q when all 10 type-2 quasars are considered: it has a modal value of 0.37 
and the region 0.23 < q < 0.59 encompasses 68% of the area. 

We summarise the derived quasar fraction as g = 0.60to;o9 for type-2 quasars showing 
narrow-lines only, and g = 0.37^*^0 for a ll type-2 quasars with spectroscopic redshifts. The 



'receding-torus' model of iSimpson I (120051 . marlced 'S05'), was derived from type-1 and type-2 
AGN showing emission lines and should only be compared to the solid blue line of Figure 4. 
His model predicts a value of g = 0.63 at high AGN luminosities, in good agreement with 
our results. 

We stress that the uncertainties are still large. For example, if instead we derive the 
quasar fraction at all z > 1.7, we estimate 6.7~^22 type-1 quasars, and the two quasar 
fractions as g = 0.45!to;i6 (using only type-2 quasars with narrow lines and 2;spec > 1.7) 
and g = 0.32+°;}^ (using all type-2 quasars with Zspec > 1-7). However, for the range 
1.7 < z < 2.6, where our spectroscopic completeness is highest, our best estimate of the 
radio-quiet quasar fraction is g = 0.37^q;^4 or, equivalently, 63^22% of radio-quiet quasars 
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are obscured. 



5. Comparison with other high-redshift samples 



Searches for PAHs in high-redshift (unobscured) quasars have proven difficult: iMaiolino et al, 

( 120071 ) observe a sample of 25 optically bright z ^ 2 — 3 quasars, find no detection even in 



their stacked spectrum, and derive an upper limit on the luminosity of the 7.7 fim PAH, 
L7.7 < 4.6 X 10^'^ Lq. The sources comprising the sample of Maiolino et al. (2007) have 
bolometric luminosities significantly larger than those of our own sample, making the PAHs 
difficult to see over the strong continuum. However, this is a good comparison sample as it 
is unbiased in terms of the star-formation rates in the host galaxies. 



Lutz et al.l (120071 ) do find a detection in the "Cloverleaf" quasar, with L7 7 = 7.6 x 10 



10 



Lq, although this is a biased case since it was observed because of it being hyper- luminous 
in the far- infrared (even after correction for gravitational magnification). 

In Figure 5, the PAH strengths and estimated bolometric luminosities of subsamples A, 
B and C are compared to those of the unobsc ured quasars mentione d above and the sample 
of submillimetre-selected galaxies (SMGs) of IValiante et al.l (120071 ). The SMGs represent 
the most powerful starbursts known, an d have their bolomet ric luminosities approximated 
to their total infrared luminosities (from IValiante et al.ll2007l ). Detection of the PAH is not 
an artefact of any selection criteria other than SMGs being galaxies with vigorous star- 
formation. 



Our sample of obscured quasars is around the break in the LF (ICroom et al.ll2004l ). so 
these sources represent the bulk of the energy density (or accreted mass density) of SMBHs. 
However, they are radio-intermediate (with Li.4ghz ~ 10^^ WHz~^sr~^), and therefore are 
more rare in space density than genuinely radio-quiet obscured quasars. In addition, this 
radio selection could possibly affect the observed properties, something that must be kept in 
mind. If the quasars were obscured solely by the torus of the unified schemes, the luminosity 
of the PAHs from star-formation in the quasar host galaxies should be similar for unobscured 
and obscured quasars. 

In our sample, emission from PAHs is detected in 3 out of 18 objects observed (subsam- 
ple A). For this subsample, the inferred values of L7 7 are similar to those of the SMGs, and 
only slightly lower than the value of the Cloverleaf. The limits inferred for non-detections 
(both for subsample C and the unobscured quasars of Maiolino et al., 2007) are similar to 
the values of the detections, so little can be said about any differences in PAH properties 
between samples. 
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The obscured quasars with detections of PAHs, have similar values of L7.7 to the SMGs, 
which are powerful starbursts. The actual conversion from L7 7 to a star-formation rate 
(SFR) is more uncertain, since for higher luminosity starbursts L? 7 does not increase linearl y 
with the far-infrared luminosity (IHaas et al.l I2OOII : iKlaas et al.l 1200 ll : iDale fc Heloul |2002| ) . 
Although it is difficult to estimate the star-formation rates (SFRs) from I/7.7, the values of 
L7.7 estimated for obscured quasars are similar to those of SMGs (see Figure 5), which have 
typical SFRs ^1000 Mq yr~^ (of all stars, integrated over a Salpeter initial mass function) 
estimated from their rest-frame FIR luminosities. These values of L7 7 are, in addition , 
larger than the limits inferred for lower-redshift X- ray absorbed quasars (Sturm et al.ll2006l ). 
or the values of L7.7 found in nearby "PG" quasars (ISchweitzer et al.ll2006l ). which have lower 
far-infrared luminosities and thus inferred SFRs (< 100 M0 yr~^) than the SMGs. 



6. Conclusions 

Two main conclusions can be extracted from this work. Firstly, we have found that the 
mid-infrared spectra are continuum dominated, and 12 out of 18 spectra show a deep silicate 
absorption feature (with r9.7~l-2). We have confirmed spectroscopically the existence of 
sources with mid-infrared spectra characterstic of heavily obscured quasars, but that do 
not show rest-frame ultraviolet emission lines (observed in the optical). The number of 
these sources is similar to the number of obscured quasars that do show ultraviolet narrow 
emission lines. Using the new mid-infrared redshifts, we confirm that the population of 
radio-intermediate obscured (type-2) quasars outnumbers the unobscured population (type- 
1), and that 631:1% of ^ ~ 2 radio-intermediate quasars are obscured. This had already 



been suggested in iMartfnez-Sansigre et al.l (120051 ). and in this paper we confirm the result 



spectroscopically: amongst radio-intermediate sources, the population of obscured quasars 
is responsible for the bulk of the accretion onto supermassive black holes. 

Secondly, we have found that obscured quasars are sometimes hosted by galaxies under- 
going vigorous star formation, with detectable PAHs comparable in strength to those found 
in SMGs. The accompanying SFR is expected to be very large (^100 M0 yr~^), so that it 
can form a significant fraction of an galaxy in ~ 10''"^ yr. The dust and gas fuelling 
this star formation are likely to be responsible for at least part of the obscuration in some 
of the sources, with the accreting black hole deeply e mbedded. Th ese properties are those 



expected from an obscured phase of black hole growth (jFabiarull999l ). and could explain why 



in some obscured quasars the rest-frame ultraviolet lines are not detectable. 
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A. New optical spectra 



For AMS18, the redshifts derived fro m optical and mid-infrared spe ctroscopy disagree 
completely (compare Tabled] to Table 1 of iMartinez-Sansigre et al.ll2006bl ). and AMS15 had 
not been observed in optical spectroscopy. AMS03 was found to have an interesting source 
nearby: in imaging data, the source was point-like in K, but galaxy-like in R, suggesting a 
reddened quasar. For these reasons the 3 objects were observed using the ISIS spectrograph 
at the William Herschel Telescope in July 2007. 



In the spectrum of AMS18 reported in lMartfnez-Sansigre et al.l (j2006bl ). three very weak 
lines were visible, at 4680, 5652 and 7527 A, yet our new spectrum shows no hint of these 
three lines. Thus, it is most l ikely that they were not real and the redshift reported in 
Martfnez-Sansigre et al.l (l2006bl ) is incorrect. Thus, we favour the redshift of z = 1.6 from 
mid- infrared spectroscopy, and this value is quoted in Table 1. 

In AMS15 we find no evidence for any optical lines, in either the blue or red spectra. 
For AMS03, the slit was placed at a different position angle (31 degrees) to go through the 
reddened q uasar. In the ID spectrurn of AM S03, we see the same double-peaked emission as 
reported in lMartfnez-Sansigre et al.l (l2006bl . Figure 6). The optical spectrum of the candidate 
reddened quasar looks like that of an elliptical galaxy at z = 0.7, and no obvious emission 
lines are visible. 

In addition, for AMS05, a new optical spectrum and imaging (using both narrow and 
broad bands) suggest the Ly a line had been missidentified as C IV, and the object is actually 
at z = 2.850 rather than z = 2.017 (D. Smith, private communication). 

Comparing optical and mid-infrared spectral properties, and radio properties, we see 
that none of the sources with fiat radio spectra show narrow lines in their optical spectra 
(AMS07, AMS15, AMS18 and AMS19, see Table [1]). This is important as, from unified 
schemes, one expects fiat-spectrum radio sources to have broad hues (except in the case of 
blazars, where the beamed continuum can outshine the broad-line region). Thus, if these 
sources are obscured, the obscuration must occur due to dust on a larger scale (not due to 
the torus) and so the narrow lines, as well as the broad lines, are expected to be obscured. 
Alternatively, the dust might be on a small scale similar to that of the torus, but with a 
covering fraction ~ 1, so that the narrow-line region cannot be excited by photons from the 
central engine. 
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Table 1: Measured and derived properties of our sample of obscured quasars. 



Name 


RA 




Dec 








Optical^ 


Radio^ 


A^l'li 




Sub— 




(J2000) 










Spectrum 


Spectrum 






sample^ 


AMSOl 


17 


13 


11.17 


+59 


55 


51.5 


- 


B 


SS 


10 


10 


N.F. 


AMS02 


17 


13 


15.88 


+60 


02 


34.2 


1.8 


B 


ss 


10 


10 


B 


AMS03 


17 


13 


40.19 


+59 


27 


45.8 


2.698 


NL 


SS 








N.O. 


AMS04 


17 


13 


40.62 


+59 


49 


17.1 


1.782 


NL 


ss 


10 


10 


B 


AMS05 


17 


13 


42.77 


+59 


39 


20.2 


2.850^ 


NL 


GP 


3 


3 


C 


AMS06 


17 


13 


43.91 


+59 


57 


14.6 


1.8 


B 


SS 


6 


6 


B 


AMS07 


17 


14 


02.25 


+59 


48 


28.8 


- 


B 


FS 


10 


10 


N.F. 


AMS08 


17 


14 


29.67 


+59 


32 


33.5 


1.979 


NL 


- 


9 


9 


B 


AMS09 


17 


14 


34.87 


+58 


56 


46.4 


2.1 


B 


SS 


8 


6 


B 


AMSIO 


17 


16 


20.08 


+59 


40 


26.5 


- 


B 


ss 


10 


10 


N.F. 


AMSll 


17 


18 


21.33 


+59 


40 


27.1 


1.6 


B 


SS 


10 


10 


A 


A 1\ /TO -1 O 

AJVlblz 


17 


18 


22.65 


+59 


01 


54.3 


2.1Q" 


ATT 

NL 


SS 








N.O. 


AMS13 


17 


18 


44.40 


+59 


20 


00.8 


1.974 


NL 


SS 


3 


2 


C 


AMS14 


17 


18 


45.47 


+58 


51 


22.5 


1.794 


NL 


SS 


3 


2 


c 


AMS15 


17 


18 


56.93 


+59 


03 


25.0 


2.1 


B 


FS 


10 


10 


B 


AMS16 


17 


19 


42.07 


+58 


47 


08.9 


4.169 


NL 


GP 








N.O. 


AMS17 


17 


20 


45.17 


+58 


52 


21.3 


3.137 


NL 


SS 


8 


6 


A 


AMS18 


17 


20 


46.32 


+60 


02 


29.6 


1.6-^ 


B 


FS 


10 


10 


B 


AMS19 


17 


20 


48.00 


+59 


43 


20.7 


2.3 


B 


FS 


5 


5 


C 


AMS20 


17 


20 


59.10 


+59 


17 


50.5 




B 


GP 


10 


10 


N.F. 


AMS21 


17 


21 


20.09 


+59 


03 


48.6 


1.8 


B 


SS 


10 


10 


A 



"The redshifts are all sp e ctrosco pic. Redshifts with 3 decimal places are from optical spectroscopy, from 



Martinez-Sansigre et al.l (|2006ar ) except AMS05 and AMS18 as discussed in Appendix lA} while those with 



only 1 decimal place are from mid-infrared spectroscopy. 

^Summary of optical spectroscopy properties. B stands for blank spectrum, NL for narrow lines. 
■^Summary of radio spectral properties from iMartfnez-Sansig re et al.l |2006bh . SS stands for steep spectrum, 

FS for flat spectrum and GP for gigahertz-peaked source. 
■^Nlli and Nll2 are the number of cycles used in each of the LLl and LL2 orders respectively. All cycles have 

a duration of 120 s. AMS03, AMS12 and AMS16 where not observed with IRS due to their high redshift. 

Tor AMSOl, AMS07, AMSIO and AMS20, it was not possible to identify securely any features and a redshift, 
so they are labelled as having no features (N.F.). Obviously, AMS03, AMS12 and AMS16 cannot be classified 
into any subsample, since they were not observed (N.O.). For the rest of the sources, the letter of which 
subsample each source belongs is quoted. 

^As discussed in Appendix |3 these two redshifts have been revised. The new redshift for AMS05 is from op- 
tical spectroscopy (D. Smith, private communication), while for AMS18 it is from mid-infrared spectroscopy. 



-17- 



Table 2: Estimated PAH strength and significance 



Subsample 


Lhoi 


T 0, 


Signif.'' 


T a 


Signif.'' 














A 


0.9 


2.2±1.8 


3.4 


4.6±1.8 


6.4 


B 


0.8 


<5.4 


<3 


<5.4 


<3 


C 


2.8 


<6.5 


<3 


<6.5 


<3 



"The estimated uncertainties include both the uncertainty in removing the continuum and the noise of the 
residual. 

''For detections we quote the significance of the lines above the noise of the residual, the uncertainty in 
removing the continuum is not included. For non-detections, we quote a 3 cr limit. 
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Fig. 1. — Individual mid-infrared spectra of all 18 sources observed. As a guide to the 
eye, a the root-mean square noise estimate is plotted as a dashed-dotted histogram. It was 
obtained by extracting a spectrum from a blank region in a 2D spectrum, and scaling it to 
the appropriate number of cycles. For sources with a spectroscopic redshift, this is quoted in 
the top right corner. The features used to derive a redshift are also indicated at the observed 
wavelength. 
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Fig. 2. — Mean mid-infrared spectra of the subsamples. Overlayed for reference are high- 
quality low-rcsohition mid-infrared spectra of four local ultra-luminous infrared galaxies: 
Arp 220 (a starburst-dominated ULIRG), Mrk 273 and Mrk 231 (both AGN-starburst com- 
posites) and IRAS 12514+1027 (with an AGN dominated spectrum). These are normalised 
to go through the data around 7 /xm. Subsample A (top): mean spectrum of 3 sources with 
both the 7.7 and 6.2 /xm PAHs visible. The wavelength of these two PAHs, together with the 
8.6 /im one, are marked with lines. Subsample B (middle): this includes 7 objects with an 
excess around 8 /im, due to the lack of 6.2 /im PAH, these are not included in subsample A. 
The mean spectrum shows the excess emission to be very broad, and the 6.2 /xm PAH is still 
not detected, suggesting this is due to the minimum in optical depth around 8 fim and not 
the 7.7 fim. PAH. Subsample C (bottom): this consists of 4 sources with clear continuum, 
and no hint of an excess around 7.7 or 8 /im. The mean spectrum is well described by that 
of IRAS 12514. 
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Fig. 3. — Residuals after subtracting the template of IRAS 12514 to the mean spectra of 
subsamples A (top), B (middle) and C (bottom). For reference, the mid-infrared spectrum 
of M82 is overlayed, which shows the expected positions of the 6.2, 7.7 and 8.6 fim PAHs. 
Subsample A has two regions of excess flux, compared to either template, these are at the 
exact wavelengths of the 6.2 and 7.7 /im PAHs, and have the correct width. Subsample B 
shows a broad excess centred at 8 fim, and no hint of the 6.2 fim PAH. Although the latter 
could be lost in the noise, the excess at 8 fim is probably best explained by continuum. 
Therefore, we do not consider subsample B to have a detection of the 7.7 /im PAH. Sub- 
sample C shows no residual, so the template provides a very good description. The region 
of 5.5-8.5 fim has lower noise than the region 8.5-11.5 /im, and of the two regions only the 
former is used to estimate the significance of the detections in Table [21 
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Fig. 4. — Posterior probability distribution function of the quasar fraction, g, in the redshift 
range 1.7 < z < 2.6, given our data and given our modelled number of type-1 quasars. 
The solid blue line shows the distribution if we only use type-2 quasars showing lines in the 
optical spectra (corresponding to rest-frame ultraviolet), while the dashed red line shows 
the case where spectroscopic redshifts from Spitzer IRS have been used to complement the 
optical redshifts. Note that only spectroscopically-confirmed type-2 quasars are used. The 
small number of sources used to derive the solid blue curve leads to a large error, so that 
the probability of g > 1 is not zero. For the dashed red line, the number of sources used is 
large enough to avoid this problem. The vertical line labelled 'S05' marks the prediction of 
the Simpson (2005) receding torus model, which should only be compared to the solid blue 
line. 
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Fig. 5. — Luminosity of the 7.7 fim PAH line L7 7 vs bolometric luminosity Lboi (both in solar 
luminosities) for samples of z ~ 2 quasars and submillimetre-selected galaxies (SMGs) with 
comparable mid-infrared spectroscopy. Upper limits are 3 cr. Our sample of z ~ 2 obscured 
quasars is represented by two points: the mean value for detections (Figure 2 panel A) and 
the upper limit from the non-detections ( Figur e 2 panels B and C). The unobseu red quasars 



are fr om the sample o f iMaiolino et al 



fl2007h . 



the "Cloverleaf" quasar is from iLutz et al. 



(120071 ). the SMGs from IValiante et al.l (120071 ). T he vertical dashed l ine represents the break 



of the unob seured quasar luna inosity function (LF lCroom et al.ll2004j ). converted to Lboi using 
the SED of lElvis et al.l (119941 ). This should only be compared to quasars (both obscured and 
unobseured) but not to SMGs. 



